On the blow up of supercritical solutions of the Nordheim 

equation for bosons. 
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Abstract. In this paper we prove, on the one hand, that the solutions of the isotropic, spatially 
homogeneous Nordheim equation for bosons, for which the initial data have bounded values of 
particle density and energy density in the portion of the phase space where Bose-Einstein conden- 
sation takes place, cannot be classical solutions globally in time. Therefore, blow-up in finite time 
in the norm takes place for this class of initial data. On the other hand, we prove the existence 
of weak, initially bounded solutions of the isotropic, spatially homogeneous Nordheim equation for 
(-H I bosons, that condense in finite time (i.e. for which a Dirac mass at the origin appears in finite, 

positive time). 

Key words. Nordheim equation, dissipation of entropy, finite time blow-up, Bose Einstein con- 
densation. 

1 Introduction. 

In this paper we continue the study of the blow-up properties of the solutions of the homogeneous 
Nordheim equation initiated in [3]. This system describes the dynamics of a dilute homogeneous 
quantum gas of bosons. We will denote as F{t,p) the distribution of particles in the space of 
moments. The evolution of F is given by the following system of equations (cf. 11 ): 



dtFi^ f I f q (F) Md'p2d''p^d''pi , pi € 
Jm? Jk3 Jr3 



, t>0 (1.1) 



i^i(0,p) =Fo(p) , pieW (1.2) 



|2 



q{F) = q:,{F)+q2{F) , e = Ml (I.3) 
M = M [pi,P2\P3.,Pi) = 5 {pi +P2-PZ - Pa) 5 (ei + £2 - £3 - £4) (1-4) 



(73 (F) = F3F4 (Fi + F2) - F1F2 (F3 + F4) (1.5) 
q2 {F) = F3F4 - F1F2 (1.6) 

^ Dcpartamento de Matematicas. Universidad del Pais Vasco UPV/EHU. Apartado 644. E-48080 Bilbao, Spain. 
E-mail: miguel.escobedo@ehu.es 

2 Basque Center for Applied Mathematics (BCAM), Alameda de Mazarredo 14, E-48009 Bilbao, Spain. 

^Institute of Applied Mathematics, University of Bonn, Endenicher AUee 60, 53115 Bonn, Germany. E-mail: 
velazquez@iam.uni-bonn.de 



1 



where we use the notation Fj — F{t,pj), j e R'^. Notice that, since we consider homogeneous 
distributions, the density F measures the number of particles for unit of volume, i.e., the number 
of particles with moment in the cube [p,p + d'^p] in a volume V would be given by F (p) V(Pp. 

In the case of isotropic distributions the system (jl.ip - (jl.6p can be rewritten in a simpler form. 
The isotropy of the solutions means that F (t,p) — F {t, TZp) for any TZ e SO (3) , p G R"^, t > 0. 
Then, there exists a function / = f {^,t) where e is as in (|1.3p such that f{t,e) — F{t,p) which 
solves: 

Q_2 roo pQO 

dth^-^J^ i, '^^^)^*3de4 (1.7) 

where: 

min{^,^, 7^,7^1} 
W= '-^ ^h=^ — - > e2 = e3 + e4-ei 1-8 

and q{-) is as in (jl.3|) with £2 = £3 + 64 — 61. More details about this computation as well as 
additional information about the physics of the Nordheim equation can be found in [3j . 

One of the main interests of the study of (jl.ip - (jl.6p or its isotropic counterpart (|1.7p is the 
relation of these systems with the dynamical formation of Bose-Einstein condensates (cf. [5] , [6] , 
[12], [l3]). It can be argued, on physical grounds, that the steady states of p.ip . (jl.3p - (ll.6p are 
the Bose-Einstein distributions: 

Fbe (p) = moS [p ~ po) + -r- ^ (1.9) 



exp(^^£^+a) -1 



where mo > 0, /3 G (0,oo] , < a < 00 and a ■ toq = Oj Po G K^. The precise sense in which 
the distributions Fbe are steady states requires further clarification, because the right-hand side 
of (II. ip is not defined in general if F contains Dirac measures. The main reason to consider the 
distributions Fbe as the steady states of (II. IL (ll.3p - p.6p is that these distributions maximize the 
entropy of the system of bosons under consideration for a given value of the momentum, energy 
and number of particles of the system for unit of volume (cf. (0])). 

One of the most peculiar features of the steady states Fbe in (|1.9P is the possible presence 
of a macroscopic fraction of particles at the value p = po- This feature is known as Bose-Einstein 
condensation. We will assume in the following that the total momentum of the system is zero. This 
can be always assumed choosing a suitable incrtial system. In this case the thermodynamics of the 
system of bosons can be described by means of two quantities, namely the number of particles by 
unit of volume M and the energy for unit of volume E : 



M = jFBE{p)d^P , E = j FBE{p)^^<i'p 



It is possible to associate unique values of toqj /? and a, with a ■ rriQ = to any pair of val- 
ues M > 0, -E > 0. We can split the set of all the values for the particle density and energy 
{{M,E) : Af > 0, E > 0} in two different phases, namely the set of values for which the corre- 
sponding value of mo in (jl.9p is positive, and those values of M, N for which mo = 0. We will say 
that the first class of states has a Bose-Einstein condensate, while the second class of states does 
not have a condensate. It is not hard to see (cf. Subsection 13.11) that the set of values {M,N) for 
which there is a nontrivial condensate is characterized by: 

M>^^(^^El (1.10) 
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where C (•) is the Riemann zeta function. 

The theory of Bose-Einstein condensates described above is meaningful for particle systems 
described by the equilibrium distributions Fbe in (|1.9p . On the other hand, the equation (jl.ll) - 
p.6p can be used to describe the evolution of a large class particle distributions Fq satisfying some 
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general condtions which ensure that the particle and energy densities remain constant for arbitrary 
times. It has been suggested in the physical literature (cf. [S], [B], [12], [13]) that the onset of 
a macroscopic fraction of particles at the lowest energy level (i.e. a condensate) is related to the 
blow-up of the solutions of (jl.lL (jl.3D - (jl.6p . More precisely, the numerical simulations in [B], [12) ') 
show the existence of isotropic solutions of (jl.7|) . with initial data / (0, e) = /o (e), that become 
unbounded in finite time. This has been proved rigorously in [3j where the authors have shown that 
there exist a large class of bounded initial distributions /o for which the corresponding solutions 
of (|1.7p blow-up in finite time. To wit, a blow-up criterium has been obtained which states that 
for any initial particle distribution with a large concentration of particles with small energy, the 
corresponding solution of (|1.7p becomes unbounded in finite time. 

We now precise the concept of subcritical and supercritical data. Given a distribution of 
particles /o we can evaluate the corresponding number of particles and energy respectively by 



We will say that /o is supercritical if the inequality (jl.lOp holds. If, on the contrary, we have 



supercritical distribution, the corresponding stationary distribution having the same amount of 
particles and energy would have a nontrivial condensate. 

We will prove in this paper that any solution of Nordheim equation, initially bounded and 
supercritical, must necessarily become unbounded in finite time. 

An interesting fact that is worth noticing is that the quadratic terms in (|l.ip . (cf. (ILGD l that 
are in principle lower order terms, play a crucial role in the onset of blow-up. The reason is that 
they produce a transfer of mass from energies of order one to small energies. 

We now sketch the main idea in the proof of the blow up result of this paper. The key 
ingredient is the blow-up criterium in P]. There are two main conditions needed to apply such 
blow-up criterium. First we need to have / > > in some average sense for small values of e. 
On the other hand, we need to have an amount of mass of order p^* , for some d^, > 0, in a interval 
[0, p] , with p small. In this paper we prove that every bounded solution of the Nordheim equation 
satisfies these conditions for sufficiently long times. To this end we use mainly two arguments. We 
first prove that the quadratic terms of the Nordheim equation transport an amount of mass towards 
small values of the energy, and that this transport cannot be balanced by the corresponding loss 
terms. As a consequence we obtain an amount of mass of order vR^ for some j/ > in any interval 
[0, -R] for sufficiently long times. The second argument that we use is the existence of an increasing 
entropy for the Nordheim equation. The corresponding formula for the dissipation of the entropy 
can be used to prove that, at least along some sequences of time, the distribution / approaches 
one stationary state of the Nordheim equation. The total energy and number of particles of the 
admissible stationary states obtained corresponds to the supercritical regime of the Bose-Einstein 
condensation. This implies that there exists a positive amount of mass in a small interval [0,/9] 
with yO > small for times sufficiently long. The criterium in implies then blow-up in finite time 
for the solutions of the Nordheim equation. 

We also remark that the blow-up criterium is local. Therefore, it is not difficult to obtain initial 
data /o ix) yielding blow-up in finite time for subcritical initial data. 

A different concept of solution for (II. 7p has been introduced in , [TU] where a class of classical 
solutions of (|1.7p which behave as / (e, i) ^ a (t) e~ e as e — ^ has been obtained. Such solutions 
have a decreasing number of particles. 

A theory of weak solutions for this equation has been developed in the papers [7], [SI. In 
those papers a definition of measured valued solutions has been given and it was proved that 
such solutions converge asymptotically to one of the isotropic equilibria distributions in (jl.9p . 
Uniqueness of these solutions for a given initial data is unknown. In the last part of this paper 
we prove the existence of weak solutions of (11.71) , whose initial data are bounded and such that a 



means of: 





subcritical. Notice that given a 
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Dirac measure at the origin appears in finite and positive time. This may be seen as the occurrence 
of Bose Einstein condensation. This finite time condensation for weak solutions of p.7p . (II. 8p is 
proved with the same method as in the proof of the blow up result in [3] , namely a contradiction 
argument based on the assumption that the initial data is a measure such that the measure of the 
origin is zero. 

The plan of this paper is the following. In Section 2 we recall the concept of solution which 
we will use as well as the Local Well-Posedness Theorem proved in 0] . Section 3 states the Main 
Blow-up result proved in this paper. Section 4 summarizes some results proved in [3] which will 
be used in the Proof of the Main result of this paper. Section 5 contains one of the key ingredients 
of this paper, namely, a lower estimate for the amount of mass contained in regions e e [0, R] with 
R small. This estimate is due to the effect of the quadratic, Boltzmann terms in (jl.7p . Section 
6 contains a rigorous proof of the entropy dissipation formula for the solutions considered in this 
paper. Section 7 contains a technical Lemma that estimates the maximum amount of particles 
contained in regions e > R, with R small, in terms of the total energy of the solution. Section 
8 proves that for the so-called supercritical data, the distribution of particles g {t, e) contains a 
positive amount of mass in the region where e is small, at least for some subsequences of time. 
Section 9 contains the proof of the main blow up theorem of this paper. Section 10 explains how 
to apply the result in 3 in order to obtain blow-up for the solutions of (ll.lOp for some subcritical 
distributions of particles. In Section 11 we prove the existence of weak solutions with bounded 
initial data such that a Dirac mass at the origin develops in finite and positive time. 



2 Well-posedness in the class of functions bounded as jy^^ 
with 7 > 3. 

A local existence result ? 

We recall the well-posedness result obtained in [3]. We first precise the concept of solution of 
pTT)) . pT^ that we wih use. 

2.1 Definition of mild solutions. 

Given 7 € R we will denote as L°° (1 + e)'') the space of functions such that: 

11/11 L-(B+-(i-h£)^) = sup {(1 + e)"^ / (e)} < 00 

^ ^ ' ' e>0 

Notice that L°° (R+; (1 -t- e)^) is a Banach space with the norm || •|lioo(]j+.(i_|_g)T) • Given T2 > 
Ti > 0, we define L^^ ([Ti, T2) ; L°° (R+; (1 + e)'^)) as the set of functions satisfying: 

SUp||/(i,-)llL~(R+;(l+£)-') < 00 

for any compact K C [Ti,T2) . Notice that these spaces are not Banach spaces. We also define the 
space ([Ti, ; L°° (R+; (1 + e)''')) which is the Banach space of functions such that: 

ll/llL°°([Ti,T2];L°°(R+;(l-l-e)T)) = II / (^i II L°= (R+;(l+e)^) < 

te [Ti ,T2] 

Definition 2.1 Suppose that 7 > 3, T2 > Ti > 0. We will say that a function 

feL^,{[T,,T2);L°° (R+;(l + er)) 
is a mild solution of |_?. 7| ), U.8\) if it satisfies: 



/(t,ei) =/o(ei)vI/(t,ei)-f ^ 



fVr^ r r/3/4(l + /i + /2)W^&3de4ds 
Jo ^(s.ei) Jo Jo 
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a.e. t g [Ti,T2), where: 

a{t,ei) = -^J^ f2{l + f3 + h)Wde3dei , * (i, ei) = exp a (s, ei) j . (2.1) 

Remark 2.2 Since 7 > 3 there exists a constant C > such that, for any t e [Ti,T2): 
hh (1 + /i + /2) t^de3de4 < C (1 + ll/ILo.(R+.(i+,)-,)) 



00 /-OO 



T/ie term a {t, ei) is bounded by C + ll/!lL~(R+-(i+e)^) j Io° Io° f2Wde3de4. By the definition of 
W , and using €2 as one of the integration variables, we estimate f2W de^de^^ as 

/o°° /2 (-y/e2 + V^) Wde2 < C ||/|litxj(R+.(i^_£)7) ■ Therefore, all the terms in \2.1]) are well defined 
for T,<t<T2tffe L^^ ([Ti, Ta) ; '(R+; (1 + e)^)) . 

2.2 Local well posedness. 

The following well-posedness result has been proved in [3]. 

Theorem 2.3 Suppose that fo £ L°° (K+; (1 + e)'') with 7 > 3. There exists T > depending only 
on ll/o (')llL°°(K+-(i+e)"') '^"''^ ^ unique mild solution of ( [j.7| ), (013) in the sense of Definition [KT[ 
/eL^J[0,T);L- (M+;(l + e)^)). 
The obtained solution f satisfies: 



47rV2^"/o(e)e"'de = 4^V2^°° /(*,£)£"' A, t e (0, T) , e | i ^| . 



(2.2) 



The function f is in the space W^'^ ((0, T) ; L°° (IR-+)) and it satisfies JJ. 7] ) a.e. e E M+ for any 
t G (OjTinax) , with initial datum f (-Tt) — /o (■) ■ Moreover, f can be extended as a mild solution 
of i-?. 7| ), il.8\) to a maximal time interval (0,T,„ax) with < Tmax < 00. If Tyanx < 00 we have: 

lim sup ||./(i,-)llLoo(R+) = 00 (2.3) 

3 Statement of the main result. 

3.1 Subcritical and supercritical data. Phase portrait in terms of the 
energy and the density of particles. 

For isotropic distributions the steady states Fbe in (|1.9I) have po — and reduce to: 

Fbe [p] = Fbe {p; a, (5, mo) = m^S [p) H — j-^ (3.1) 

exp + -1 

with mo > 0, /3 e (0, 00], < a < 00 and a ■ niQ = 0. (3-2) 

It is customary in the physical literature to denote a as —(3fi, where < is a quantity 
with units of energy termed as chemical potential. We will use the notation (j3.2l) in order to use 
nonnegative quantities in the arguments. If a = and mo = the resulting distributions are the 
Planck distributions. 

The following result shows that at equilibrium there exists a one-to-one relation between the 
values of the particle density and the energy and the values of the chemical potential, temperature 
and density of particles in the condensate state. 
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Proposition 3.1 Given Af > 0, -B > there exists a unique steady state 
Fbe (p; o:, P, mo) in the family i3.S\) such that: 



FBEip;a,l3,mo)(fp = M 
We have mo — iff 



[ FBE{p]a,!3,mo)^-d?p^E 



c(l) (^V 



M < — . (3.3) 

(c(i))- V3; 

where C (•) is the Riemann function. 

Proof. This result is well known and it can be found in classical texts of Statistical Physics (cf. 
for instance [Ij, [^j). We sketch here the main arguments. Suppose that mp = 0. Then, since a > 
we have (cf. ((L9)) ): 

M < / j-^^ = Stt - / — — = (3.4) 



exp(^)-l \PJ Jo 

a 2 

where we have used spherical coordinates, the change of variables well as Taylor's series 

to expand (1 — e^^) ^ . A similar computation yields: 

1 f M^d^p /2V f°° xidx 

E = - r— TTT = TT ' 



2 _ 1 \PJ Jo e-~l 



Then, if mn = we obtain M < ^^'^ , (4f ) ^ £'5 . On the other hand, if mn > 0, we can argue 

as in the derivation of p. 41) and using the fact that a = we obtain M > (^^^^C(|)- Combining 
this formula with (13.41) it then follows that M > — -i (^) ^ E's , whence the Proposition follows. 

■ 

Definition 3.2 We will say that a couple of values {M,E), such that M > and E > 0, is in 

the subcritical region if iS. 'J^) holds. On the contrary, if M > — 3 {-r)^ E's we will say that 

(Cd))^ 

(M, E) is in the supercritical region. 

3.2 Blow-up Theorem for supercritical data. 

The main result of this paper is the following. 

Theorem 3.3 Suppose that fo G L°° (M+; (1 + e)^) with 7 > 3. Let us denote as M, E the 
numbers: 

l*QO 1*00 

An / fo (e) V2^d€ ^ M , in fo (e) V2^de = E (3.6) 
Jo Jo 
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Let us denote as f e L^J[0, T^ax) ; (R+; (1 + e)'')) the mild solution of [T^ , ifTT^j m 
Theorem \2.3\ where Tmax is the maximal existence time. Suppose that: 

(c(i))' V3; 

Then: 

Tm^K < oo (3.8) 

Remark 3.4 Notice that M, E are the particle density and energy density associated to the dis- 
tribution fo- The numerical factors are due to the change of variables to spherical coordinates and 
to the fact that we use e instead of p as an integration variable. 



Notice that Theorems 12.31 and Theorem 13.31 imply that for supercritical values of (M, E) , f 
blows-up in finite time: 



Corollary 3.5 Suppose that the conditions of Theorem \3.3\ hold. Then: 

lim sup |!/(i, •)IIl°=(k+) = °o ' T„iax < oo 



4 Summary of some results in [3]. 

We now recall some results which has been obtained in [3] which will be used in the Proof of the 
results of this paper. 



4.1 Blow-up criterium. 

One of the main ingredients in the Proof of Theorem 13.31 is the following criterium for blow-up of 
mild solutions which has been proved in [J. 

Theorem 4.1 Let M > 0, E > 0, > 0, > 3. There exist p = p{M,E,v) > 0, K* = 
K* (A/, E, v) > 0, To = To (A/, E) and a numerical constant > independent on M, E, v such 
that for any fo G (IR^; (1 + e)'') satisfying 

47r / fo (e) V2^de = M , Att [ fo (e) V2?d€ = E (4.1) 

/ fo{e)V~ede>i^R^^ for Q < R < p , (' fo{e) ^ede> K* {pf' (4.2) 
Ja Jo 

there exists a unique mild solution f € Lf^^ ([OjTmax) ',L°° (M+; (1 + e)'^)) defined for a maximal 
existence time T^ax < Tq. The solution f satisfies: 

lim sup ||/(-,OIIl-(r+) = oo 

Remark 4.2 Theorem \4.1\ has been formulated in f^jj with the first inequality in {^-/^ replaced by 
/o(e) > V on (0, R) for R sufficiently small. A careful examination of the Proof of Theorem \4-.l\ in 
JSjl shows that only the inequality in j^.iip is used in the proof of the result. 
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4.2 Transfer of mass formula. 

We will use the following result to estimate the rate of transfer of particles between different 
regions. This is a combination of several results obtained in [3]. 

We recall that the density of particles for unit of energy is given by (cf. [3]): 

5(t,e)=47rN/2^/(t,e) (4.3) 

Then: 

Proposition 4.3 Suppose that 7 > 3 and f E ([0, T^ax) ; L°° (1 + e)'')) is a mild solution 
of I-?. Tp , d j.<^|) in the sense of Definition \2.1l Let g be as in i4.3^ . Suppose furthermore that 
ip £ Cq ([0, T) X [0, 00)) Then, the function %j:^p (t) ~ Jj^^^ g {t, e) ip (t, e) de is Lipschitz continuous 
in t E [0,T] and the following identity holds: 



dt(^Jg(pde^^ Jgdt(pde+-^ j j j '^"^1111 ^^ (ei, £2, £3) deide2de3~ 



]_ f f f 919293'^, 
^ f f f 5132$, 



2 2 Jr+ Jr+ 7r+ veleiei 
+ ^ I I I ^i^H^&ideadea , a.e. t e [0,T] (4.4) 

i+ Jr+ Jr+ V^1^2 



with: 



$ = min I y/e^, y/e^, ^ {ei + £2 - £3) + } (4.5) 
'H^ = if> (£3) + (ei + £2 - £3) - (ei) - V (£2) (4.6) 
^/^ (ei, £2, £3) = g X] (<^'T(i),edr(2),e(T(3)) * (ecr(i),edr(2);e(T(3)) 

where is the group of permutations o/{l,2, 3}. Moreover, for any convex function ip we have: 

ev(ei,e2,e3) >0 , (£1, £2, £3) e R+ (4.7) 
G^{ei,e2,£3) =Gv{ea{i),<^cr{2),<^cr{3)) for any a eS^ (4.8) 

Proof. It is just a consequence of Lemma 3.13 and Theorem 4.2 in [3]. ■ 

5 A lower estimate for the mass in [0, R] with R small. 

Proposition 5.1 Suppose that fo G L°° (R+; (1 + £)''') with 7 > 3 and let M, E he as in h3.b]) . Let 
f, Imax the mild solution of Jj. 7| ), H.S}) and the maximal existence time whose existence is proved 
in Theorem\2^ respectively. Then, there exist Tq ^ Tq {E, M), p = p (E, M) and K = K{E, M) 
such that j/Tinax > Tq {E,M), we have: 



I g (£, t) de > KRi for any < R < p 
Jo 



for any t > Tq {E, M) 



Proof. We will assume that p < 1. Given < i? < 1, we define the following family of test 
functions: 



RJ 



+ 



where (s)_|_ — max{s, 0} . Using Proposition 14. 31 we obtain: 

dt( j gVRde] > 5 / / / ^^^H^^deide^dea , a.e. t e [0,T] 
\J J ^ Jr+ Jr+ Jr+ v^i'^2 
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where we use the fact that ipa is convex, and then (|4.7|) holds. Using (I4.6P as well as the fact that 
if (ei + £2 — 63) > we obtain for a.e. t e [0, T] : 



dt (^J g(pB.de^ > h - h 



where 



h = - 
2 



h = - 



2 Jb+ JR+ iR+ \/£l<^2 



<i5i?. (£3) dtxdtidt-i 

{Vr (ei) + </5fl (£2)] deide2de3 



(5.1) 

(5.2) 
(5.3) 



We estimate using the definition of $, as well as the symmetry with respect to the permu- 
tations of ei, £2, as follows: 



31.92* 



R 



< TT / 5l(p_R (£1) d£l / g2d€2 

Jo Jm+ 



(fiR (£l)d£ld£2d£3 

min l^eT, ^} (fi + £2) 



/£i£2 



where we use that $ < min {^/^i, and we compute the integral in £3. Notice that 

min{^^ie,+e,) ^ ^ ^^^^^ ^.^^^ £l < i? < 1 : 



92 



R+ 



rin{^, (£i + £2) 



/£i£2 



rf£2 < 2 / g2max{l,^/€^}d£2 

<2 I 32 (1 + £2) rf£2 < 2 (M + £;) 



/2 < 27r {M + E) / .9^fld£ 
Jo 



(5.4) 



We now estimate Ii assuming that gip^de < Using the fact that J gde — M we have: 



gde > 



M 



(5.5) 



Indeed, suppose that gde < then J^^ gde > whence, since ipR > \ for < £ < 



we have gfRde > > ^ against the assumption, whence (|5.5p follows. 



Using the definition of $ (cf. (14. Sp ) we obtain: 



where: 



( ' fl ) 



+ Jr+ V^1<^2 



9i92deide2 



(5.6) 



J{Vi,V2)= min (771 + 772 - 773) + ! (1 - 773)+ d% (5.7) 

Notice that, if 771 > i, 772 > 5 we have J (771, 772) > cq > where cq is a number independent 
on R. Then: 



CO /'OO 



/ ('11 

V i? ' r) 



/£l£2 



gig2deide2 



> 



CO /"CO 



£ JH \le\e2 



3152 , , 7rcoi?2 
-.de\de2 — — - — 



d£ 



(5.8) 
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Due to the boundedness of the energy we have: 



Then, ii L > ^ , gde < ^, whence, usmg ^ 



^ , M 
i 16 



Therefore (jS.Sp yields: 




512i 

We have then obtained that there exists a = a {E, M) such that, 

h>aRi. (5.9) 
Combining (15.11) . (|5.4p and (|5.9p we obtain: 

9t sV'iide) > - 27r [M + E) gifade if ^ gipnde < ^ (5.10) 

Choosing R < p < ^liM^lM^j jt folfows that, for gipude < ^^"^^^^^^ we would have also 
Jq^ g'^nd^ < X ^^"^ inequality (j5.10p holds. This implies that J gipude increases exponen- 

3 

tially until reaching the value ^^"^^^^^ or larger in a time which depends only on M and E. If 

Jo 9fRd^ > the inequality J gipjide > ^^"^^^^^ is automatically satisfied, due to our choice of 

p. On the other hand, if 4^^"^^^) — lo^ 9^Rd£ < x value of /J* g^PRde remains in that region 
due to ([SHI)) . Therefore, there exists Tq = Tq (M, £;) such that if Tq (M, £;) < T^^^: 



aR2 



47r (M + E) 
whence the result follows. ■ 



< J g^Rdt < y gde, for t > To (M, E) 



6 Dissipation of entropy formula. 

We now formulate the formula for the dissipation of the entropy. The entropy associated to the 
bosonic, isotropic quantum Boltzmann equation (jl.7D . p.8p is given by: 

S[f]^ [ [(l + /)log(l + /)-/log(/)]V^de (6.1) 

We also define, for any given / the dissipation of entropy by means of: 

D[f]= [ [ [ (l + /l)(l + /2)(l+/3)(l + /4)[Ql,2"Q3,4]x 

JR+ JR+ JR+ 

X [log (Qi,2) - log {Q3a)] $deide2de3 (6.2) 
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where: 

£4 = £1+ £2-63 , $ = min (ei)_^, ^ (ea)^, -y/ (£3) + , (^4)+} 



We then have: 

Proposition 6.1 Suppose that fo £ (K+; (1 + e)''') with 7 > 3, and /o ^ 0. Suppose that 
f G ([0, Tmax) ; L°° (R+; (1 + e)'')) is the mild solution of [U^ , IT^) mi/i iniiza/ dotom fo as 
in Theorem \2.3\ Then, there exists C — C {E,M) such that: 

\S[f{-,t)]\<C{E,M) ,0<t<r^ax (6.3) 

For any Ti and T2 such that 2Tq < Ti < T2 < T-^^x, where Tq as in proposition i5.1\) . we have: 

0<i^[/(-,t)]<C(Ti,r2)<oo , Ti<t<T2 

for some constant C {Ti,T2) . Moreover, the following identity holds: 

S [/] (T2) - S [/] (Ti) = r D [f (., t)] dt (6.4) 

Proof. As a first step we prove (|6.3p . To this end we use the following inequalities: 

[(l + /)log(l + /)-/log(/)] <C/ if / >! 

[(l + /)log(l + /)-/log(/)] <C/|log(/)| if 0</<l 

We then can estimate 5* [/] for any / G (K+; (1 + e)'') , 7 > 3 as follows. We have the 
inequalities: 

S[f]^ f [■■ ■] V~^de + f [...]^ede+ [ [■■ ■] V~ede 

J{f>i} "'{/<l: |l0g(/)|<V^} J{f<l, |l0g(/)|>V^} 

<C [ fVede + C [ feUe + C [ 77 (77 |log (/)|) 

J J J{/<cxp(-e)} ^ ^ 

<C{M + E) + C J cxp (-| j de = C (M, E) 
whence (16.31) follows. 

We now need a lower estimate for f {e,t) , i > in order to ensure the convergence of the 
integrals in (|6.2p . To this end we use methods inspired in the ideas used by Carleman to prove 
the 7?— Theorem for the classical Boltzmann equation (cf. [5], Section 5). We claim that for any 
w > and < Ti <T2, there exists Cti,T2,ui such that: 

/(£,i) > CTi,T.,c.exp (-£i+") ,Ti<t<T2 (6.5) 

where the constant Cti,T2,uj depends also on /q. 

In order to obtain (|6.5p . we first remark that the mild solutions of (II. 7p . (II. 8p satisfy (cf. [3], 
Lemma 3.7): 

Q_2 /-oo ^00 

at/l + ^Myil/i > /3/4 (1 + /l + /2) T^d£3d£4 

o_2 pco pco 

> ^ / / fsfiWde^dei (6.6) 
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for a.e. t e [Ti,T2\ and where Cq > 0. 

Since fo^O there exists an interval with the form (i?i , ORi ) with i?o > and 6 > 1 perhaps 
close to one, such that /^^ ^ fode > 0. Moreover, due to the continuity of the mild solutions in time 
in the weak topology in space it follows that, assuming that t is small we have: 

f {t,e)de>mi > , 0<t<t (6.7) 

-Ri 

Let us choose i] > satisfying (2 — ry) < 2 , (2 — 77) > 6* > 1. We can then derive a lower 
estimate for if ei G ((2 — 77) (2 — i]) 9Ri) . Combining (I6.6P and (|6.7p we obtain the estimate: 

dtfi+T:M^fi>K^,e(^J^ ' fde^ , ei e {eR^ {2 - rj) OR^) (6.8) 

Iterating this argument we obtain the chain of inequalities: 

dtfi + nM^h > K,^e (^j ^ " fd^ , ei e ((?i?„, 0i?„+i) (6.9) 

with: 

i?„+i = (2 - 77) i?„ . (6.10) 

Let us consider the sequence of positive times: t„ = t — t„ with t„ = with < 6 < 1, 
independent of n. Using the fact that / > we obtain, using Duhamel's formula: 

fiit,ei) > Kn^e exp {-t: M y/e^ {s - tn)) f{s,e)dej ds 

for a.e. €1 G {eR,„0Rn+i) (6.11) 

If we assume that b is small, depending only on i, we would obtain that t„ > 0, for all n > 1. Using 
(|6.11|) we then obtain: 

fit,ei)>K^^gaJ^ (^J^ "/(s,e)dej ds , ei G (0i?„, 0i?„+i) (6.12) 

for some a > independent on n. Applying (j6.7p and (j6.12l) with tt, = 1 in the interval (i?2, dR2) C 
{9Ri, OR2) , we obtain: 

fi{t,ei)>K,,^ga{t-t2){mif , ei € (i?2, ^^^2) , a.e. t e [t2,t\ (6.13) 
Integrating this inequality in {R2,0R2) we obtain: 

/ f {t,e)de>m2^ K^.ea{e -l)Ri{t~t2){mif , a.e.te[t2,f\ (6.14) 

Due to the definition of the sequence {tn} we have: 



(^3 - t2) = (t-T3- ia) = (t2 - T3) = — E= I 1 -^^ ) = ^^"^ ''1 [i^J- 

Using (|6.14p with t — ts it then follows that: 

— R 
m2 > K^,e^a,b—7==ml > Cr,fi,a,b.Rt'm\ 
V R2 
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Iterating the argument it then foUows that: 



i+i > C'7,,e,a,f),i?,i (2 - r/) 2 m„ , n = l,2,3 



Moreover, since the term (2 — ri)~ increases to infinity we have that, after a finite number of 
steps the factor in front of becomes larger than one. We obtain lower estimates for to„q in 
terms of the first term. Then: 

TO„+i >m1 , n>no 



Iterating we then obtain: 
Since m„Q is a fixed number, in general smaller than one we obtain: 



rrin > {runaf ° , n>no 



m„ > exp /3i2" ' 

with /?! > independent on n. Since _R„ is of order (2 — r/)" this implies an estimate of the form: 

m„ > Cexp(-^i?i+") , n = l,2,3,... 

for some /3>0, C>0,aj>0 independent on n, where moreover ui is such that 1 + lu > 

Notice that, since 77 > may be taken as small as we wish by taking 9 sufficiently close to 1, the 

value of u! may be as small as we need. 

Using now (|6.12l) as well as the fact that the union of the intervals 6'i?„) , n > 1 cover the 
whole interval (i?i,cxi), we obtain: 

/(t,e)>CTi,r.,c.exp(-ei+") , e > i?i , < Ti < t < T2 (6.15) 

where the value of w may have been changed from respect to its previous value in order to eliminate 
the constant (3 from the exponential, but still remaining as small as we need. The constant Cti,T2,uj 
depends on /q. 

In order to obtain a lower bound for small values of e G (0, Ri) we use again (|6.6D . By (|6.15p . 
for any ei > we have: 



Z poo poo 



g 2 roo .ei 



2Ri Jo 

El 



> / V^fsdes. (6.16) 



Using now Proposition (|5.1I) we have, for all ei E (0, p{E, M)): 

Jo 

From (I6.16P and (|6.17p we deduce 

dtfi + TrM^h > Csi, for t > Tq, 
and therefore, after integration in time from to t for any > Tq: 

fit, si) > CVT, ( 1 - e--^'^v^(*-^-)) , for t > T,. 



(6.17) 
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If moreover {t — T^ ) < 1: 

fit,ei)>Ceit>CiTi,T2)ei. 
This shows that for some positive constant Cti,T2,uj' 

3T 

/(i,£i) > Cti,t.,^£i, eie(0,p), -Y <Ti<t<T2. 

where p is as in Proposition (jS.ip . Arguing as in the proof of formula (6.15) above, we deduce after 
a finite number of iterations that for some positive constant Cti.T2,u'- 



.f{t,ei)>CT^.T,,u£u eie(0,i?i), 2To<Ti<t<T2, 

where the constant Cti,T2,(^ rnight change from hne to hne. 
We deduce then from (j6.15D and (j6.18p : 

f{t,ei)>CT,,T2,u:£iexp{~e^+-^), £i > 0, To{E, M) < Ti < t < T2. 

Using (|6.19p we obtain: 

|log(/)| <Cti,t. .^(1 + log £ + 6^+") , e>0 , 0<ri<t<T2. 



(6.18) 



(6.19) 



(6.20) 



We now use (j6.20p to show that the different integral terms appearing in the formula of the 
dissipation of the entropy (cf. (j6.2p ) are finite. Due to the boundedness of / we just need to show 
that the following integrals are finite: 



/1/2 
/1/2 
73/4 
73/4 



log 
log 
log 
log 



/l 



1 + /1 
h 

1 + /3 
/l 

1 + /1 
h 

1 + /3 



log 
log 
log 
log 



/2 



1 + /2 

/4 
I + /4 

/2 
I + /2 

/4 
I + /4 



$c?eide2(ie3 
^deide2de3 
^deide2de3 
^deide2de3 



where £4 — ei + £2 — £3- Replacing the variable £4 by the variable £1 in the last two integrals, 
and relabelling the number of the resulting integration variables (namely £2, £3, £4 to £1, £2, £3) we 
reduce the estimate of the last two integrals to the first two ones. Using now the symmetry of the 
variables £1, £2 we are left only with the three different terms: 



h = 

/2- 
/3 = 



/1/2 
/1/2 
/1/2 



log 
log 
log 



/l 



1 + 
h 

1 + /3 
h 

1 + /4 



^deide2de3 
^deide2de3, 

$d£ld£2(i£3 



Using (|6.20p and the boundedness of / we obtain: 

Il<C f I I /1/2 (1 + £l+") $d£irf£2d£3 



c 



{ei>e2} 



[• • •] $d£id£2d£3 + C 



^deide2de3 



{ei<e2} 



14 



We now use that $ < in the first integral and $ < ^JT\ in the second. Then, since 
£3 < ei + £2 

h<C I f /1/2 (1 + el+^) e,^deide2+ 

J J{ei>e2} 

+ C f f f /1/2 (1 + e^y/^deide^ 

Since, at is has been indicated before, the value of ut in (j6.20|) may be chosen as small as we 
need, we will assume in the following that w < (7 — 3) . Then, using that / < (^i^^y ■ 

h<C [ [ {1 + e,)-'' {1 + e^y {1 + el+^) ei^de,de2+ 

J J{ei>e2} 

+ C [ [ [ {1 + e,)-'< {1 + {I + el+^) e2V^deide2 

J J J {ei<C2} 

< c 

In order to estimate I2 and I3 we use a symmetrization argument that yields: 

h<C f f f (1 + ei)-^ (1 + £2)"^ (1 + + log £3) $&i&2de3 

J J J{£i>e2} 

<C [ [ [ (1 + (1 + (1 + ey + logei) eiV^deide2 

J J J{ei>e2} 

< C 



h<C (1 + ei)-^ (1 + 62^ (1 + el+" + log £4) ^deide2de3 

J J J{£i>e2} 

<C f f (1 + ei)-^ (1 + £2)-^ (1 + + log £1) ^de^de2 

J J{ei>e2} 

< c 

where in the estimate of I3 we have used that £4 < 2ei. 

In order to conclude the proof of (j6.4p we need to use a symmetrization argument. We will use 
that (cf. 3 , Theorem 3.4): 

= ^ / / W [(1 + /i) (1 + /2) fsU - (1 + /3) (1 + U) /1/2] deader , a.e. t e [Ti,T2] 
V 2 J J D(ti) 

Differentiating (|6.1I) and using the conservation of mass, we obtain the following: 

a(5|/l) = ///.og(i^)*x 

X [(1 + /i) (1 + /2) /3/4 - (1 + /3) (1 + k) /1/2] ^£1*3*4 

log 



X [Q3,4 - (9i,2] $deide3de4 



We now claim that for any function / e L°° (M+ : (1 + e)''^) satisfying (|6.15p the following 
identity holds: 

Ji - J2 (6.21) 
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where: 

Ji= I I / log 



X [(93,4 - Qi,2] '^deidesde^ 



J^^\j j j I [log (Q3,4) - log (Q1.2)] [Q3.4 - Ql,2] (1 + /l) (1 + /2) X 

X (1 + /a) (1 + /4) $deide2&3 

In order to prove (jOTI) suppose first that / G L°° (K+ : (1 + e)^) f] C(R+) satisfies (|615)) . 
Then: 

./l = jy y y y[l0g(Q3,4)-l0g(Ql.2)] [Q3.4-Qi,2](1 + /i)(1 + /2)x 

X (1 + /3) (1 + /4) $(5 (ei + 62 - £3 - £4) deide2de3de4 

= J2 

whence (|6.2ip holds for continuous functions. For arbitrary functions / G L°° {M.^ : (1 + e)'') 
satisfying (16.15^ we can obtain (|6.2ip approximating / by means of a sequence of continuous 
function /„ converging to / at almost every t e [Ti,r2] . ■ 

7 Reformulation of the criticality condition: A technical 
Lemma. 



We now prove an auxiliary result that reformulates the conditions in Proposition 13.11 in a form 
that only depends on the values of the equilibrium distributions for values of e > i? > 0, with R 
small. We define a class of auxiliary functions: 

/,(e;a,/3) = \ , e > -a , /3 > 0, aeR (7.1) 

exp [p (e + a)j - I 

The following result holds: 

Proposition 7.1 Given E^, > 0, and S > there exist Rq {E^,, S) > 0, Lq {E^,, 6) > such that, if 
fs is one of the functions in i7.1\ ), < R < Rq {E^, S) , L > Lq {E^, 6), a > and 

E = 4tt [ /,(e;a,/3)V2^de , E < E, 

JR 

then: ^ 

47r / f^{e-a,[i)V2ede<^^^(^\ ET. +5 (7.2) 

Proof. We will assume in all the following that i? < 1, L > 2. Suppose that E < E^,. We define 
a family of functions function Pe,r,l • i ^ by rneans of the relation: 

E = 4n [ f, (e; a, Pe,r,l {a)) V2?de 

JR 

Since the functions fs (e; a, /3) are strictly decreasing with respect to /? and, for any given 
a > — we have: 







47r lim 


f 




Ir 






4tt lim 


f 




Ir 
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it follows that the function Z^^^ij^i is well defined. We now define the functions Af£;_ij;L : [~"fjOo) 
R+ as: 



Me,r,l (a) = 47r / fs (e; a, Pe,r,l (a)) V 2ede 



We now claim that 



dME,R,L {a) 
da 



< 



for a e [— -f 7O0) ■ Indeed, differentiating Me,r,l we obtain '^^'^^^q^^""' = ^^[^| , where: 



Ai (a) = Svr / / dy 
Jr Jr. 



g^i(x+a)g/3(y+Q) 



(e^*(^+") - 1)^ (e/3fe+") - 1)^ ' 



^xy^ (x + a) — \/yx^ {x + a) 



(7.3) 



A2 (a) 



da; 



> 



'i?^ (e/3(^+") - 1)^ 

Symmetrizing the variables x, y in Ai (a) we obtain: 

Ai (a) = ~8tt / dx / ^ ^ ^ ^ < 

Jr Jr (e^(^+") - 1)' (6/3^+") - 1)' 



whence (|7.3p follows. Therefore: 



M 



E,R,L 



(a) < 



E,R,L 



R 



where 13e,r,l ) satisfies: 



/2ede 



(7.4) 



E^A^I^f.^e^-l^E^R., (-f 



Then, since it! < 1, L>2 

c2 



R 



^de 



4-77 



/2e3de 



It then follows from the fact that E < E^, that there exists /3* — /3* (£'*) > such that: 

R 



Pe,r,l ( ) >/3* 



(7.5) 



(7.6) 



(7.7) 



Using (|7.7I) it follows that, for any Eq > 0, there exists Lq — Lo(-E'*i£o) > such that, if 
L>Lo: 

4^y^ fs[e;^^,PE,R,L[-^y^V2^de<eo (7.8) 



We define the functions: 



$fl (/3) = 47r /, (^e; - 1, \/2;^de 



(7.9) 



We define E by means of: 



E = E + A^j^ f,^e;-^,PE^R,L(^-^ 



'2e^de 
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Using (TTS]) we obtain E < E < E+Eq. Notice that due to ([73]) we have $_r {Pe.r.l (-f )) = E. 
We now claim that there exists £i > depending only on i?* , 5 such that E < ei and a > — 



we have Me,r,l (a) < 6. Indeed, (j7.6p implies that Pemx 



can be made arbitrarily large if 



£1 is small. Then (j7.1|) and (|7.4p implies that Me,r,l (a) can be made smaU due to Lebesgue's 
dominated convergence Theorem. Then (j7.2l) would follow in this case. We remark that £i is 
independent on R. 

We will assume then in the following that E > ei. We claim that there exists /3* — (3* (i?*, £i) , 
such that Pe,r,l {—^) < P* ■ Indeed, we can estimate i>_R (/3) as: 



^R (/3) < 4^ 



(e^^ - 1) 



-de < Stt 







and the right-hand side of this formula converges to zero as /3 oo due to Lebesgue's dominated 
convergence Theorem. Then, f3E,R,L (-f ) < /3* = /?* (E^^Si) ii E > ei > 0. We now claim that 
the functions (/3) converge uniformly in the interval [/?*,/?*] as i? to the function: 

V2^ 



Jo Jq 



{el^'^ - 1) 



de 



(7.10) 



Moreover, we have also uniform convergence in the interval /3*] of the derivatives {(3) to 
% as i? ^ 0. 

Indeed, this just follows from the inequalities: 



e + 



'de + C / y/^d< 



\^R (/?)-<j>o m<c 



l*'i?(/?)-1>o(/?)l<C 



Since $0 (/?) is strictly negative for /3 e [/?*,/?*] we then obtain that, if E > £1, (3e,r,l (— "f) 
satisfies: 



-/^•"^de + C / Vide 



E,R,L 



HE) 



< Cen 



(7.11) 



if R is sufficiently small, where (3 [E) is the unique solution of: 

$0 {E)) = E 

Similar computations yield the uniform convergence in the interval j3 G /?*] of the functions 
{(3) defined by means of: 



^b,{(3)=Att f,(^e;-^,(3^V2'ede 



to: 



^-0 W) = 47r 



2e 



(e'^^-l) 



de 



if > £1. Using (|7.4p we can write: 



Me. 



R,L 



Pe,r,l 



+ *o 13e,r. 



IV 



(7.12) 



(7.13) 
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Due to (I7.11|) and the uniform convergence of the functions 'i'n to '^q stated above, it foUows 
that the term between brackets in (j7.13p can be made arbitrarily small if R is small enough. On 
the other hand, using again (17. lip we can make 4'o {(3e,r,l arbitrarily close to {P {E}) 

if Eq is small. Therefore, if R is small we obtain: 

Me,rx (^-|) < *o (/9 (E)) + 5 (7.14) 
Using the definitions of $o and ^Pq in (j7.10D . (|7.12p as well as Proposition 13. II we have: 
*. (^)) < (I) ' ,^,))» ^ (f ) * iE)i 

Combining this estimate with (|7.14p the result follows. ■ 



8 Mass concentration in the region of small energies for a 
sequence of times tn — > +00 



We will prove now that if the solutions of (11.71) are globally bounded and (j3.7p holds, the corre- 
sponding functions g (t,-) would have a significant amount of mass in the regions where e is small. 
The main result of this Section is the following. 



Proposition 8.1 Suppose that /o, / are as in Theorem \3.3\ Let us assume that Tj^ax — oo. Then, 
there exists > and p > 0, both of them depending only on M, E such that, for any < R < p 
there exists a sequence {tn} with f„ — > oo as n oo such that: 



/ g{t„,e)de = ATT y^ef {1^, e) de > (8.1) 
Jo Jo 



for any n. 

In order to prove Proposition 18. II we need several Lemmas. We begin deriving an estimate for 
the number of particles with large energy. 



Lemma 8.2 Suppose that fo and f are as in Theorem \3.3[ Then, for any Eq > 0, there exists 
L = L {E, £q) such that 



/ g {t, e) de = 47r / %/2^/ (t, e) de > M 
Jo Jo 



forte [0,r,„ax]. 

Proof. It is just a consequence from the conservation of energy E as well as the inequality: 



/oo 1 /"^ 

9 {t, e)de<j-J^ g (i, e) ede 



E 

T 



Choosing L > — the result follows. ■ 
We define the following auxiliary function: 

Q{t,e)- 

Thcn: 



l + /(t,e) 
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Notice that: 

0<g(t,e)<l , Q{t,e)<f{t,e) , e > 0, te[0,T„,ax] (8.3) 

Then: 

/OO 
V2^f{t,e)de<eo , 

with L ^ L {E, £o) as in Lemma [8.21 Moreover 

/"OO /"OO 

47r / V2^Q {t, e) de < M , An V2^Q {t, e)de<E (8.4) 
Jo Jo 

We define also: 

«'(s) = slog(l + s) (8.5) 
We will use the following concept of weak convergence. 

Definition 8.3 We say that a sequence {Qn} C L°° (R^) converges weakly to Q*, and we will 
write in this case Qn Q* iff-' 

lim / Qnfde — / Q^(pde (8.6) 



n— >oo 



+ 



for any test function G Co [0, oo) . 

Remark 8.4 // the sequence of functions {Qn} satisfies < Qn < 1, we can apply a density 
argument to show that i8. 6\) holds for any ip £ (0, oo) . 

We have the following result: 

Lemma 8.5 Suppose that fo, f are as in Provosition \8.1i There exists a sequence {t„} , i„ — > oo 
as n —> oo such that: 

III f {tn,e,) f {tn, e,) ^ f ^ ^ l^^'^l - l) <^de,de,de, -> (8.7) 

as n ^ 00, where £4 = ci + £2 ^ £3. Moreover: 

Q (tn,-) Q* {■) as 00 (8.8) 

where £ L°° (M+) and < (e) < 1 , e > 0. 
Proof. Notice that ((^ yields: 

^ 1^' ^ L L L * - 

Due to (ESI) we have |S'[/](r2)| + |5[/](Ti)| < 2C{E,M). Then, since T^ax = 00, and 
l-S* [/] (72)1 < 00, we can use ((O)) to obtain: 



/"OO 

/ D[f{-,t)]dt<^ 



Therefore, there exists a sequence , i„ — > 00 such that D [f (-jtn)] — > as 71 — > 00. Using 
then (|8.9p we obtain (|8.7p . Using (|8.3p and classical compactness results for measures in the weak 
topology, we can then extract a subsequence of {tn} (which will be denoted in the same way) for 
which ([STg]) holds. ■ 
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Lemma 8.6 Suppose that fo, f are as in Lemma \ 8.5l Let us assume that Q^, — in i8.8\) for 
0<i?i<e<i?2<oo. Then: 



R2 



Ri 



as n ^ 00 



(8.10) 



-Ri 



Proof. Choosing the test function ip such that (/3(e) = 1 for e E [i?i,i?2] and <p (e) =0 for 
e ^ [Ri, R2] (cf. Remark l8.4p . we obtain, using the nonnegativity of Q, that for n large enough: 



rR2 

/ Q [tm e) \/erfe — > as n — 00 

jRi 



In order to prove the second formula in (j8.10p we define the set: 



{R1.R2) = <! i?i < e < i?2 : Q [tn, e) > 2 



Then: 



/(in,ei)/(in,e2)* 



Q (in, £3) Q (in, £4) 



> 



U„(RuR2) JU„{Ri.,R2) "'{Q(t,>,C3)<i} 



Q (in, El) Q (in, £2) 

/(in,ei)/(in,e2) x 



1 ^deide2de3 



^ ^j, f Q{tn,e3)Q{tn,e4) \ 

VQ (in,ei)(y (in,e2) / 

Using that Q (i,i, 64) < 1, Q {tn,ej) > \ if G Wn (-Ri, ^2) for some j = 1, 2, as well as the fact 
that the function ^ (s) is decreasing for s e (—1, —5) we obtain: 



, e2) VI, f Q(/»^^3)Q(in,e4) _ \ ^^^^^^^^^^ 



> 



^^'"''^^^^'"''^^^ VQ(in,.l)Q(in,e2) 

/ / /(in,ei)/(in,e2) X 

UARuR2) Ju„{Ri,R2) J {Qit„,<^3)<i}n{^<e3<^} 

X * (4Q (t„, 63) - 1) $de3deide2 



> 



U„{Ri,R2) 



f{tn,e)de] X 



{Q(t.,E3)<4}n{4L<.3<^} 
\ 2 r 



*(4Q (i„,e3)-l)cie3 



fitn,e) de 



U,ARi,R2) 



Due to the first formula in (|8.10p we have that 



Q{tn,e3)> <e3< 2 



Ri 



y{Q(t..«3)<|}n{-^<e3<^} 



3i?i 



dez 



as n — > c». Then: 



lim 

n— foo 



{Q(t„,.3)<i}n{i^<.3<^} y{i^<.3<^} 
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Therefore, using (|8.7p : 

lim / f{tn,e)de^O 
We then have, taking into account ()8.2p and the first hmit in (j8.10|) : 



hm / /(t„,e)v^de= lim / f {tn,e) ^/ede+ 

f Hm / /(t„,e)de< 2 hm / g e) T^de = 0. 



Lemma 8.7 Suppose that fo,f are as in Lemma \8.5[ Let us assume that ^ for e > Ri in 
I18.8\) . Suppose that < i?i < 1. Then: 



\Q (<n, £3) Q {tn, £4) - Q {tn, El) Q {tn, £2)! deide2de3 ^ as n ^ 00 (8.11) 

/fli JRi 

Proof. Using (18.31) wc obtain: 

+ (in,£l)(y (in,£2) / 

<l I I f {tn, £l) / £2) ^ i!"' ^ i!"' - l) $d£ld£2rf£3 



Q (in,£i) Q (i«,£2) 



Then: 



> 



/ {tn, £l) / £2) * ( Q(;"'^3)Qfa,£4) _ ^\ ^^^^^^^^^^ 

rrf^nu \n(, ^,,.f Q{tn,e3)Q{tn,e,) A ^, , , 

/ / / Q{tn,ei)Q{tn,e2)^ [-r— -— - 1 \ <i>deide2de3 

J Br, Jri \Q{tn,ei)Q{t„,e2) J 



^^rrf^nu .,..f Q{tn,e3)Q{tn,ei) A,,, 

- \hr / L Q{tn,£i)Q{tn,e2)'i> \-pT7-, 77777 r - I ] deidezdea 

V ^ Jri \Q{tn,ei)Q{tn,e2) J 

We define a convex function ^ {s) as: 

§ (s) = * (s) if s > 0, and * (s) = * (-s) 

Notice that 4' (s) is bounded as C |s| for large \s\ . Using Jensen's inequality we obtain: 

2 



/ f°° \ Ft - 



'Ri 

Iri Iri I^l \^ *3 (^»' - Q (^n, £1) Q {tn, £2)! dtide2de^ 



<l I I f {tn, £l) / £2) ^ i!"' ^ i!"' - 1 ) 1'&1&2&3 



,Q(in,£i)(3(i«,£2) 

Taking into account (|8.7p as well as the fact that ^ (s) vanishes only for s = we obtain (|8.1ip . 
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Lemma 8.8 Given < i?i < i?2 < oo, suppose that fo as in the statement of Proposition 
\8.1l let {tn} as in Lemma \8.5l Suppose that ^ m the interval [Ri,R2] . Then (e) = 
exp (— /?* (e + a)) for e € -R2] with > and a > —Ri. Moreover, there exists a subsequence 
of {tn} which will be denoted luith the same indexes such that: 

Q{tn,-) ^ Q*{-) tnL^{Ri,R2) as n 00 . (8-12) 

Proof. Let us define: 

i?* = supjr G [i?i,i?2] : J Q*de = o| (8.13) 
Since ^ in the interval [Ri, R2] we have < i?2- We then define also: 

/•CO poo 

Qnies)^ / \Q {tn,€3)Q {tn,ei + €2 ~ €3) - Q {t„,ei)Q {tn,e2)\d€id€2 

J R JR 

where R = niax{i?i, , where 6* < 1 is very close to one, to be determined later. 

Using Lemma 15771 and more precisely (18. lip , with Ri = R, it follows that there exists a set 

V C , ^ with measure |V| — R and a subsequence of {tn} , which will be labelled with the 

same indexes, such that: 

Qn (£3) -^0 , n 00 for any £3 G V (8-14) 

The definition of R^ and the assumptions of the Lemma imply the existence of a (5 > such 
that 

/ Q,{e)de>0 (8.15) 
Jr,,+s 

We now claim that there exists G V n {i?* < e < R,^ + 6} and a subsequence of {tn}, labelled 
with the same indexes such that: 

lim Q{tn,e,) = ry > (8.16) 

Indeed, otherwise we would have lim„_j.oo Q (e,in) = for any e G V n {i?* < e < R^, + S}. 
Lebesgue's dominated convergence theorem, combined with (j8.3D would imply that (e) = 



a.e. e G 



2 , i?* + (5 



2 

Therefore: 

00 /'OO 



, but this would contradict the definition of i?* in (|8.13p . 



/ / \Q {tn,£*) Q {tn, ei + €2 - e*) - Q {tn,ei) Q {tn, £2)1 deide2 ^ as n ^ CX) 
JR JR 

Writing x — ei ~ e^,, y — €2 ~ ^* and _ff„ (x) — Q {tn, e* + x) it then follows that: 

/•oo poo 

/ / \Q{tn,e^)Hn{x + y)-Hn{x)Hn{y)\dxdy^Oa.sn-^oo (8.17) 

JR-€, JR-e, 



We define: 



Gn{x) = [ Hn{Od^ (8.18) 

Note that this mtegral is finite for any n, due to (|8.3p . Moreover, they are uniformly bounded 
by a constant C = C (i\f, Ri) < oo. Then: 

/ \Q{tn,e,)Gn{x)-Hn{x)Gn{0)\dx 

JR-e, 

/>oo />oo 

< / / \Q{tn,e,)Hn{x + y)- Hn{x)Hn{y)\dydx 

J R-t, JQ 
poo poo 

< / \Q{tn,e,)Hn{x + y)- Hn{x)Hn{y)\dydx^Q 

JR-e, JR-e, 
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where we use that e* > i?* > R. 
Defining: 

A„ (x) - Q {tn, e,) Gn (x) - H„ (x) G„ (0) (8.19) 

it then follows that: 

dx ^ as n — )• oo (8.20) 





X„ (x) 


JR-e, 





Due to we have G„ e M^^'"" (i?* - e*, oo) . Then (jSlB implies: 

A„(x) = g(i„,eOG„(x)+G„(0)G:,(:E) , a.e. x>i?-e, (8.21) 
Notice also that the weak convergence of the sequence Q (•, in) implies that G„ (0) — ?► Q* (e) de 

G„(0)' ^" ~ G„(0) 



as n oo. Due to (|8.15l) we have (e) de > 0. Let us write A„ (x) — ^" ^ , /3, — 



Due to (|8.16p . (|8.20p and the fact that lim„_^oo G„ (0) > 0, we obtain: 

/>oo 

/ |A„ (a;)|da; ^ , /3„ ^- /3* > as n ^ oo (8.22) 

JR-e, 

G; (x) + /3„G„ (x) = A„ (x) (8.23) 

Integrating (|5231): 

/•OO 

G„ {x) = A.e^'^"^ + / e-'^"(^-«)A„ (y) dy , x>R-e. (8.24) 



for suitable constants A„ S R. Due to (j8.22p the integral term in (I8.24p converges to zero, uniformly 
in the set [0,oo) . Since G„ (0) is uniformly bounded it then follows that the sequence {An} is 
bounded. Taking a new subsequence if needed, it then follows that An — > A* as n — > oo, whence: 

G„ (x) A*e"'^*^ , x>R~e^ 

Using (j8.23p we obtain that _ff„ — > A^fH^e"^'^ in {R~ for any L fixed, sufficiently 

large. Therefore 

(3(t„,e) ^- exp(-^* (e + a)) as n oo , aSR (8.25) 

in {R,L) with L large. We now consider two cases. If i?* — Ri, (|8.25D would imply (|8.12p . 
Otherwise, we choose 9 < 1 sufficiently close to one to have i?i < 9R^ < i?*. Then (|8.25p would 
imply (e) — exp (— /3* (e + a)) and this would contradict the definition of 7?* in (j8.13p . Therefore 
-R = i?* = Ri. Using that (e) < 1 for e > it follows that a > —Ri and this concludes the 
Proof. ■ 



Lemma 8.9 Suppose that fo as in the statement of Frovosition \8.1\ Let Q < Ri < 1 < R2 < oo 
and {tn} as in Lemma \8. 8\ Let us write /* (e) = cxp(ff (e+a))-i ^^^''^^ '^1 P* '^'"^ '^^ ™ Lemma \8.8[ 
Then, for a suitable subsequence of {tn} which will be labelled with the same indexes and for any 
6 > small: 

f{tn,-)^f*{-) in L\{Ri,R2)n{Q^ <1-5S}) as n ^ 00 
with as in Lemma \8.8l 

Remark 8.10 Notice that the restriction {Q* < 1 — 5(5} applies only to the limit distribution Q*, 
and not to the sequence Q (tn, •)• 

Proof. Taking a subsequence if needed we can assume that the convergence Q{tn,-) — ^ Q* (•) 
in Lemma 15^ takes place for a.e. e G (i?i,i?2). Let us denote as I the set X = (i?i,i?2) n 
{e ; Q* (e) < 1 — 56} . We estimate the L^ norm of / (t„, •) — /* (•) as follows: 

|/(i„,e) - {e)\de = Ji,„ + J2,„ 
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Ji.n= / |/(i„,e)-/46)|de (8.26) 



J2.n= / \fit^,e)-Me)\de (8.27) 
Bn^{e>0;Q{t^,e)>l-6} (8.28) 

The sequence J2,„ converges to zero as n — >■ oo due to (|8.2p as well as the fact that the function 
^ is bounded in the corresponding integration region. 



(1-Q(t„,e))(l-Q.(e)) 

To estimate Ji^„ we use Lemma 18.51 Then: 

/ / / / (tn. ei) / e,) * - <i>d63d£i ^ as n ^ oo 

Notice that $ can be estimated from below uniformly in n if, say, 

i?i < £3 < ^ (ei + £2) (8.29) 

On the other hand, since we integrate in £3 in I we need to ensure that the domain where 
(|8.29p holds has an intersection with I whose measure can be estimated from below. This can be 
seen because the values of ei, €2 must be also in the interval I. Notice that I = (e, i?2) for some 
e depending on 5, i?i . Therefore we obtain that the region of integration for £3 can be replaced by 
a set (e, ^) . Notice that this set is contained in I. Then: 

Iff f{tn,e^)f{t^,e2)^( ^fr''int'''i ~l]de3de,de2^0 as n ^ 00 
Jine„ Jin6„ J(e,f ) \Q [tn, ei) Q [tn, £2) ) 

due to the fact that we have a lower estimate for $ independent on n for £3 e (£, ^) . We now 



use the convergence of Q (t„, •) to Q* (•) in (£, which is a Corollary of Lemma Egoroff 's 
Theorem shows that there exists a set C (e, ^) with a measure arbitrarily close to | where 
Q (tn, •) — > Q* (•) uniformly. Then: 

/ / / / {tn, 61) / {tn, £2) * f n!!"''ini!"'''! - 1) ^3deide2d ^ as n ^ oo 
Jlr\B„Jlr\B„JA \^ (hi, ei) [tn, €2) / 

and assuming that is sufficiently large we would have: 

Q{tn,e3)Qit„,e^) ^ (1 - 56) < ^ _ ^^ 
Q{tn,ei)Q{tn,e2) ~ {l~Sf ~ 

if S is smaU. Then, - l) > co > independent of n, whence: 

/ f {tn,e)de\ — > as n cx) 

/ f{tn,e)de^O as n^c50 (8.30) 
JlnB„ 

On the other hand, since Q* (fi) < (1 — 56) for £ e Z, it follows that: 

|XnS„| ^0 as n ->cx3 (8.31) 
due to Lebesgue's Theorem. We can take the limit of the sequence Ji^„ as: 

Ji,n< [ fitn,e)de+ [ Me)de^O 
JinB„ JinB„. 



Co\A\ 

Then: 
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due to the fact that (e) is bounded in I and using (j8.30p . (|8.3ip . ■ 

Proof of Proposition 18.11 Given /o as in the statement of Proposition 18. 1[ we have that the 
stationary solution i^se (p; a, /3, toq) in Proposition 13.21 having the number of particles M and 
energy E satisfies mg > 0. We define m* — Due to the Lemma TS. 2 1 we can select L > such 
that: 

g{t,e)de^47r f {t , e) V2'ede < (8.32) 

for any t S [0,Tniax] • 

We now apply Lemma 1531 Then (j8.8p holds. We have now two possibilities. Suppose first that 
= in [-J, cx)) . Then Lemma [8T6l with i?i = and i?2 = L, combined with (j8.32p imply that: 

5(t„,e)de= / g{tn,e)de+ / g{t„,e)de< / g {tn,e) de + < ° 



R 



10-5 



if n is sufficiently large. Then: 



' / \ , 4mo 
g [tn, €) de > —— > m» 



whence the result follows in this case. 

Suppose now that ^ in [-f,oo) . We then set Ri — ^, R2 = L. We choose S > small 
and apply Lemma 18.91 to obtain: 



/ gitn,e)de^4TT [ (e) %/2^de (8.33) 

f g{U,e)ede^An f U (e) 

J (R.L)n{Q, <l-55} J(fi,L)n{Q.<l-5(5} 

as n — >■ CJO. (Notice that the integrations are made in the interval {R,L) in spite of the fact that 
the convergence in Lemma [8.81 is in e > In particular a > — -f ). 

We now claim that the right-hand side of (|8.33p can be made smaller than M plus some small 
error term if S and R are small. Actually taking S much smaller than we would obtain that 
{R, L) n {Q* < 1 - 56} = (i?, L) . Then: 



/ g{t„,e)de^ 4tt [ f,{e)V2ede (8.34) 

J[RX) ■J{R,L) 

f g (t„ , e) ede ^ 4n f (e) V2^de 

J(R.L) J(R.L) 



l(R,L) •J{R,L) 

We now use that J,^ ^, g {tn, e) ede < E, whence: 



47r / /, (e) V2^de < E 
J(R,L) 



■,L) 

On the other hand, using (|8.32p as well as the fact that J gde = M we obtain that: 

5(t„,e)de+ / g{tn,e)de>M 
(0,-R.) Hr^l) 

Using (I8.34P as well as Proposition 17.11 we obtain, assuming that n is large and taking 5 = 
in Proposition 17. II 

9it..^d.<^(^El^^ 



(C(|))- 
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Using now the hypothesis p.7p in Theorem l3.3l as well as the fact that M 



(C(l)) 



mo we obtain: 



(o,i?^) 



c(i) /4.y 



3 

£■5 



mo 4mo 



5 5 



(8.35) 



if R is sufficiently small, depending only on E, M. We now use Proposition 15.11 combined with 
(|8.35p and Theorem 14.11 implies that / must become unbounded in finite time. Notice that we 
apply Theorem 14. 1 1 taking as starting time Ti = The Theorem 14 . 1 1 can be applied there due to 
the invariant of the problem under translations in time. Note also that given a mild solution of 
(|1.7p . (jl.SI) defines also a mild solution in any time interval [T'i,T2] C [0,Tmax) • Indeed, we just 
need to split the integral in time and using the remaining terms as new initial data for the solution 
using the semigroup property of the exponential term. ■ 



End of the Proof of Theorem 13.31 It is just a consequence of Theorem 14.11 Proposition 15.11 
and Proposition 18.11 Notice that the blow-up criterium 14.11 can be applied starting at any time 
Ti > and not necessarily at Ti =0, due to the invariance of (|1.7I) under translations in time. ■ 

10 Blow-up for subcritical data. 

It is worth to notice that the blow-up criterium in ^ is purely local. As a consequence it is possible 
to find initial data with values of the particle density and the energy (M, E) in the subcritical region, 
but on the other hand yielding blow-up in finite time. We formulate the result as a Theorem due 
to its independent interest. 

Theorem 10.1 Given fo € (M^; (1 + e)''') with 7 > 3. Let us denote as M, E the numbers: 



Let us denote as f e L'^^{[0,Tn,^^) ; L°° {R+ ; {1 + e^)) the mild solution of (flTTp , ifTTg)) in 
Theorem \2. 8\ where T,„ax is the maximal existence time. There exist funtions fo € L°° (M^; (1 + e)'^) 
with 7 > 3 such that: 



for which T-^^y^ < 00. 

Proof. It is just a consequence of the fact that (|4.1[) . (14.21) can be obtained for initial data satisfying 
(|10.ip , as it can be seen modifying locally fo in a region e < i? in order to obtain fo (e) > v there 
as well as an amount of mass of order one in such ball, keeping constant the values of E, M by 
means of changes of fo in regions where e > 1 . ■ 

11 Finite time condensation. 

It is usual in the physical literature to relate Bose-Einstein condensation phenomena, at the level 
of the kinetic equation (|1.7p . (|1.8p with the onset of a macroscopic fraction of particles at the 
energy level e = 0. Due to this, it is natural to say that a weak solution of (|11.2p . (|11.3p develops 
a condensate if (|11.8p holds. We notice that the methods in 3_ can be used to obtain suitable 



9 End of the Proof of Theorem [33 





(10.1) 
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finite time condensation results for weak solutions of (|1.7p . p.8p . In order to use measured valued 
solution of this equation a suitable concept of weak solution must be introduced. Since we are 
interested in condensation phenomena, it is convenient to reformulate (I1.7P for the mass density g, 
instead of /. Suppose that g{t,t) = AiT^/2ef {t,e) (cf. Then, formally (|L7)) . (fTS)) becomes: 

poo poo / n \ 

dtgi = 327r3 / qi — ^ <^de^dei (11.2) 

Jo Jo V47rV2e/ 

$ = min{^/eT, V^a, \/e3, v^} > £2 = £3 + £4 - ei (11-3) 
We will denote as 7W+ (M+; 1 + e) the set of Radon measures g in R+ satisfying: 



(1 + e)g{e) de < 00 
We will use the notation g (e) in spite of the fact that g is a measure. 

Definition 11.2 We will say that g e C{[0,T);M+ (M+; (1 + e))) is a weak solution of HT^) . 
hll.S]) with initial datum go G (K^; 1 + c) i (f? for any ip E Cq ([0, T)x[0, 00)) , the following 
identity holds: 

- I 30 (e) (^, e) c^e = / / gdtipdedt + \ f f f f ^^j^^^Qipdeide2dezdt+ 
JR+ Jo Je+ 22 Jo JR+ JR+ V^i^zca 



/O JR+ 

^ r f f f 51.92$ 



2 Jo Jr+ JR+ Jr+ ^£1^2 



Q^deide2de-idt (11.4) 



where $ is as in and; 

= (£3) + (ei + £2 - £3) - 29? (ei) (11.5) 

The existence of weak solutions of (|11.2p . (|11.3p in the sense of the Definition 1 1 1 . 21 for general 
initial data fo has been obtained in [S]. 

We now have the following result: 

Theorem 11.3 Let > 0, £; > 0, > 0, 7 > 3. There exist p = p{M,E,iy) > 0, K* = 
K* {M, E, v) > 0, To — To {M, E) and a numerical constant 9^ > Q independent on M, E, v 
such that for any weak solution of U1.2\) . U1.3\} in the sense of Definition \11.2\ with go G 
7W+ (R+; (1 + e)) satisfying 

47rV2 / go{e)de = M , 47rV2 / go[e)ede = E (11.6) 

Jr+ Jr+ 

/ 9o{e)de>vRi forO<R<p , [' go i^) de > K* (pf' (11.7) 
Jo Jo 

we have: 



0<t<To 



inf / g{t,e)de>0 (11.8) 



{0} 



Proof. The proof is essentially the same as the one of Theorem 14.11 obtained in [3]. The proof 
in that paper uses a contradiction argument, where it is shown that the existence of a mild so- 
lution of (|1.7|) . (|1.8p in the sense of Definition 12.11 satisfving supQ^j^^^^ ll/(^r)lloo ^ ^ yields a 
contradiction. A detailed examination of the Proof obtained in that paper shows that the only 
arguments used in the derivation of such contradiction are the identity (jll.4p and the condition 
supQ^(<y^ g {t, e) de < 00, whence the result follows. ■ 
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Remark 11.4 It is worth noticing the specific point where the condition on supp^^j^^ J^^j g {t, e) de 
appears in the Proof of Theorem \4-l\ as well as in the Proof of Theorem \11.3[ This condition plays 
a crucial role in the crucial Measure Theory Lemma 6. 1 of paper 131. Indeed, such a Lemma holds 
only if J^^^ g {t,e) de = 0,. Using this Lemma yields the contradiction, as in J^. Actually this is 
the only place where this condition is used in the whole Proof of these results. 



Notice that the construction of the weak solutions of (|11.2p . (|11.3p docs not rule out the 
possibility of having "instantaneous condensation", i.e.: 



inf / (t, e)de> 

0<t<T' 7(0} 

for any T* > 0. It is possible to construct weak solutions of (|11.2p . (|11.3p such that 

sup / g {t, e) de = 

0<t<T. J{0} 

for some < T* < Tq, but satisfying (|11.8p . For such solutions we would have then condensation 
in a finite, but positive time, as it has been suggested in the physical literature (cf. [S], [3], [12], 
[T5]). More precisely, we have: 

Theorem 11.5 Let M > 0, E > 0, > 0, j > 3. There exist p = p{M,E,v) > 0, K* = 
K* (M, E, v) > 0, Tq — Tq (M, E) and a numerical constant 6** > independent on Af, E, v such 
that for any fo G L°° (M^; (1 + e)''^) satisfying 

47rV2 / fo (e) ^/e^^e = M , AttV2 [ fo (e) Ve^de = E (11.9) 

[ fo{e)V~ede>i^Ri for < R < p , f fo (e) V~ede > K* {pf' (11.10) 
Jo Jo 

there exists a weak solution g of hll.2]) . hll.S\) such that go (e) — 47r\/2e/o (e) ■ Moreover, there 
exists > such that the following holds: 

sup |l/(t, •)I1loc(r+) < oo , inf / g{t,e)dt>Q (H-H) 

0<t<T, ^ ' "*<t<"oJ{o} 

where g — \m\f2ef . 

Proof. We construct a weak solution of (|11.2p . (|11.3p in the sense of Definition II 1.21 as follows. We 
first use Theorem l2.3l to obtain a bounded mild solution / of p.7p . (ll.Sp in the sense of Definition 
12.11 in a time interval < t < . We then use the approach in 8, to obtain a weak solution g 
of (1112]), (ITOj) defined for < t < oo with initial datum g[T^ ^) = 47r V2e/ (T* , ) , with / as 
obtained in the previous step. We then construct a global weak solution g defined in < i < oo 
by means of g (t, •) — 47r\/2e/ (t, ) for < t < T* and g (t, •) ~ g {t, •) for t > T^. Due to Lemma 
3.13 in [J , g is a weak solution in < ^ ^ Using then the continuity of g (t, •) , in the weak 
topology for t ~ T^, it then follows that the constructed measure g G C ([0, T) ; (M+; 1 + e)) is 
a weak solution of (|11.2p . (|11.3p defined in < t < oo in the sense of Definition 111.21 It is readily 
seen by construction as well as Theorem 111.31 that g satisfies (|ll.lll) whence the result follows. ■ 
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